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We study restricted multiple three wave interaction system by the inverse scattering method. 
We develop the algebraic approach in terms of classical r-matrix and give an interpretation of the 
Poisson brackets as linear r-matrix algebra. The solutions are expressed in terms of polynomials of 
theta functions. In particular case for n — 1 in terms of Weierstrass functions. 



I. RESTRICTED MULTIPLE THREE WAVE INTERACTION SYSTEM 

Several studies have appeared recently on coupled quadratic nonlinear oscillators [H , (2| , El 

dbj 1 



ih = 0, (1) 



d£ J 2 J 
dc ■ 1 

1 ~df + u * bj + 2 ejCj = °' (2) 

du 



'— +E^' c i= ' ( 3 ) 



3=1 

where £ is the evolution coordinate and tj corresponds to the normalized wave number mismatchs. The system 
(HEJ) is introduced in Q to model the growth of a low frequency internal ocean wave by interaction with higher 
frequency surface waves and is used in [5j as a model of plasma turbulence. This system describe triads of waves 
(dj,bj,u), j = 1, ...n evolving in £ and interacting with each other trough multiple three wave interaction with 
possible applications in optics. 

II. LAX REPRESENTATION 

Let us consider coupled quadratic nonlinear oscillators 

dbj 1 

dc ■ 1 

1 ~df + u * bj + 2 ejCj = °' (5) 

*|+X>*=0, (6) 
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where £ is the evolution coordinate and tj are constants. The equations ((4][6|) can be written as Lax representation 

§ = [M,L], (7) 
d£ 



of the following linear system: 

dip 



= M(£, X)m A) L(£, A) = 0, (8) 



where L, M are 2x2 matrices and have the form 



L ^' X )~ {c(Z,\) D(£,\)J ' (9) 



m,x) = {-$ 2 $ 2 )- do) 

where 

w-«E^J. ( 12 ) 

n c 6* 

C(£,A) = a(A) [m«-,^^ ] , (13) 

where £)(£, A) = — A) and a(A) = II"=i(^ — e 0- The Lax representation yields the hyperelliptic curve K = (y, A) 

det(L(A) - = 0, (14) 



where I is the 2 x 2 unit matrix. The moduli of the curve (fT4"]) generate the integrals of motion Jo , Jj, Kj, 3 — 1, ...,n, 

^ 2 = yL 2 (£,A) + 5(£,A)C(£,A). (15) 
The curve (|15[l can be written in canonical form as 

2n+2 
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iJIfA-A^EfA), (16) 



where Aj ^ A& are branching points. From (TT5")) and explicit expressions for A), -B(£, A), C(£, A) we obtain 
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Next we develop a method which allows to construct quasi-periodic and periodic solutions of system (OS]). The method 
is based on the application of spectral theory for self-adjoint one dimensional Dirac equation with quasi- (periodic) 
finite gap potential IX = — u cf. Eqs. (14I5|) 

- m 2 j - iA 3 -*y = 0, (19) 

*-^-U**y+*Aj-*y = 0, (20) 

with spectral parameter A and eigenvalues Aj = iej/2. The equation ((TJ) is equivalently written as 

— =iuC-iu*B, A(£,A) = £A„ +1 _^)A^ (21) 



or in different form we have 



dB 

— = -tXB - 2iuA, B(Z,\)=J2 B n -j {HW , 

dC n 

— =i\C + 2iu*A, C(£, A) = £ C ^-3 > 



Aj+i£ = iuCj — iu* Bj, Aq — 1, Ai — ci, 
iBj + i = —Bj£ — 2iuAj + i, Bq = —2u, 
iCj+i = Cj£ - 2iu*A j+1 C = -2u*, 

where c\ is the constant of integration. Differentiating Eq. Q21j) and using (fT5|) we can obtain 
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Using (??) the eigenfunction 'J'i for finite-gap potential XL have the form 
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Special case of system (|4]|6|) is the three wave system 



dAi dA 2 dA 3 



The corresponding elements of Lax matrices are 



M(£,A) = 



+ -ieAi - i^A 3 ^ 



-ieA\ -i\2 

To integrate the system (f2"5|) we introduce new variable 

1 <L4 X A 3 A| 



^ = |A 2 | 2 -|A 3 



A x di Ai 



in terms of which our equations can be written as 



where 



R{\) = (Ja 2 - 1a) 2 - |^| 2 (A - M )(A - /x*) = 
-A 4 — aiA 3 + a^A 2 — a^X + a^. 



The \i variable and Aj, j = 1, . . .3 obey the equations 



ai = 0, |^i| 2 - -A = a 2 , 

-\A 1 \ 2 (n + f i*) =a 3 , -.A-/i//L4i| 2 =a 4 . 



which are related to the integrals of motion of the monomer system with = O.The equation of motion is then 

(f) 2 = -4(M 4 + ^M 2 -^) (37) 
where the system (f2"5|) conserves the dimensionless variable N and the Hamiltonian H 

N = |Ai| 2 - ^A, H = A X A 2 A* 3 + A 3 A{A 2 , (38) 



Solving Eqs. (|36|) for p, variable we obtain 



1 



^-[ff + ^j, (39) 



where 

P(u) = Av 3 - ANv 2 + N 2 v - H 2 . (40) 
We seek the solution A\ in the following form 



= vMf)exp(#«)), (41) 



where v = \Ai\ 2 = p(£+a/) + Ci, p is the Weierstrass function, and ui' is half period. Using Eq. |3"9"|) and the following 
equation 

d ^ = -2iv{ l i- i i% (42) 
derived from (|31j) and three wave equations we obtain 



whose solution can be expressed in terms of the Weierstrass elliptic functions as 



Av 6 - ANv 2 + N 2 v-H 2 , (43) 



3 

Substituting this expression in Eq. (jlTj) we obtain 



N 

v = p(£ + »') + -. (44) 



.1, = \/p(S + w >) + E exp(#(0), (45) 



where the phase is given by 

2 P '(k) V a^ + uj' + n) 
and ipo is initial constant phase. 

III. BI-HAMILTONIAN STRUCTURE 

In this paragraph we will compute r-matrix algebra of restricted multiple three wave interaction system. We note 
that in Lax representation @ we remove the function a(A), which is essential for studying Hamiltonian dynamics of 
restricted multiple three wave interaction system. 

Let as consider Lax matrix 

&{\) = a-\\)L&\)=(* ®\ . (47) 



and introduce standard Poisson bracket, {•. -}o 



If X =_ (?£d9__dldg\ _ A ( df_dg_ _ 9/5^ 



du du* du* du I ' \ dbj db* db* dbj 

3=1 \ J •* j 

\dc 3 dc* dc*dc 3 ) 
The entries of j£f satisfy to the following well known equations 

K(A),^( M )}o = {^(A),^( M )}o = {^(A),^(/i)}o - 0, 
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K(A),^( M )}o 

K(A),^(/i)}o = ^-(^(/i)-^(A)), 
{^(A),^( M )}o = ^-(V^)-^(A)), 



which may be rewritten as linear r-matrix algebra 



{J?(A),J?(/i)} = [r(A - m),-^(A) ] , (48) 



1 2 

hereof (A) = -S?(A) ® I , _£f (/x) = I (8 -2f (m) an d — m) is a classical rational r-matrix: 



r(A- M ) = y^-, n 



/1 0' 
10 
10 

\0 1. 



Remind, that two Poisson brackets {., .}o and {., .}i are compatible if every linear combination of them is still a 
Poisson bracket. The corresponding compatible Poisson tensors Pq and P\ satisfy to the following equations 

[P Q) Po] = [Po, Pi] = [A , PI = 0, (49) 

where [., .] is the Schouten bracket iGj. Remind that on a smooth finite-dimensional manifold j$ the Schouten 
bracket of two bivectors X and Y is an antisymmetric contravariant tensor of rank three and its components in local 
coordinates z m read 

{X,Yr k = - V [X mk ^— +Y mk ^— +cyde(i,j,k) 
The Poisson bracket associated with the Poisson bivector P is equal to 

{m,g(z)} = (df, Pdg) = J2P i \z) %^ d Q z Z J - (50) 

i,k 

Here df is covector with entries df /dzi and (., .) is a standard vector product. 

There are a lot of the Poisson brackets {., .}i compatible with the linear r-matrix bracket (j48|) similar to the 
quadratic Sklyanin algebra 0]. Here we consider two examples only. 

Proposition 1 // 

z— 1 i— 1 i— 1 

where hi,ei, fi are dynamical variables and e; are numerical parameters, then the following brackets are compatible 
with linear r-matrix bracket RSy 

{^(A),^(/i)}! = MA),^0*)}i=Q, 



{^(A),#Gi)}i 
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A^(/x) - /^(A)J - ^(A)^(/i), 
*f(/i) -#(A)) +«(A)ST(m), 



(52) 



- A^(A)) +2 (l - ^(A)) ^(/x), 



{^(A),^(^)}i = 2(&f(v)tf(X) -^(A)^(/x) 



Proof : It is sufficient to check the statement on an open dense subset of the linear r-matrix algebra (|48[) defined by 
the assumption that all the hi,ei,fi and ti are different. Namely, substituting rational functions (|51[) into (|48[) one 
gets canonical brackets on the direct sum of n copies of sl(2) 



{hj,ej} = ej, {hj,fj}o = -f 3 , {ej, fj} = 2hj, j = l,...,n. 



(53) 



Substituting these rational functions (|51[) into the second brackets (I5"2"j) one gets the following non-local brackets 
between generators hi , a , fi 



{hj, ej}i = (ej - ej)ej, {/ij, = -fe - e^/j, {e h fjjx = 2(e^ - e^/ij 



(54) 



ej}i = -eibj, {hi, fj}i — etCj, {e l7 f }i = -2eiOj , {fi, f j } 1 = -2hiC 3 ■+ 2f, 



Now it is easy to prove that Poisson bracket (|5"3"|) is compatible with the Poisson bracket 
Proposition 2 // 
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(55) 



where hi,ei, fi are dynamical variables and e; are numerical parameters, then the following brackets are compatible 
with linear r-matrix bracket 



{<%{\),mp)}i = w),4)h = o, 
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A^(^) - M^(A)J - pi^(A)^(/i), 
<T(m) - ^(A)) + p!^(A)^(/i) - P2^(A), 

- A.fi/(A)) + 2^1 - pi^(A)) .0^) - p 3 



(56) 



Here 



MA),%)} = -2 Pi (^(A)C(/x)-^(mMA)) +2p 2 (^(A)-^(/x)) + P3 (^(A)-^(/i)) 
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where \SC jW\ is a quotient of polynomials X and Y in variables A and /i over a field as in Q/. 



Proof: As above it is sufficient to check the statement on an open dense subset of the linear r-matrix algebra 
defined by the assumption that all the hi,ei, fi and e, are different. 

Namely, substituting rational functions l[51"j) into l[48|) one gets local brackets n — 1 copies of sl(2) 

{>',■',}■■■ {hjifjh = -fj, {ej,fj}o = %hj, j = l,...,n-l, 

and degenerate brackets 

{e„, /«} = -2h n , {h n , hi} = {h n , ej = {h n , fi} = . 



(59) 



(60) 



The leading coefficients h n is the Casimir element for these brackets. 

Substituting rational functions (|5ip into the second brackets (151)1) one gets second non-local brackets between 
generators hi, ej, fi. At i, j = 1, . . . , n — 1 these brackets looks like 



{hj,ej}i 



{hjJsh 



2 e 



(61) 



{hi,ej}i = J -, {hi,fj}i = — J -, {e l Jj\i = -. 

f-n &n C n 

At e n = 1 these brackets coincide with the previous brackets (j54"l) . 
The remaining non-zero brackets have the following form 



{fh fj}l — 



-2h 



2fiQj 



{fi, fn}i = P3(A + M = ej) /•*, {e», / n }i = -P3(A + /u = £j) ef, {e n , /„} = 

Now it is easy to prove that Poisson bracket (|59 |) - (fHU|) is compatible with the Poisson bracket ijBTj) - 
the proof. 

In the both cases we can rewrite second Poisson brackets in the following r-matrix form 

{if (A), Jf ( M )}! - [n 2 (A,/i),r(A)] - [r 2 i(A,/i),T(/x)], 



where 



f*i2(A,M) 



/a^ 

^ 



\ o o o 



X—fi 



0' 



-p 3 

,P2-piC(Ai) 0, 



(62) 

This completes 
(63) 

(64) 



and 



r 2 i(A, /i) = nri 2 (M,A)II. 
For the Lax matrix Jzf with entries ([51]) we have 

Pi = 1, Pi = 0, p 3 = 0. 



For the Lax matrix j£f with entries (|5"5f functions pk are given by 1|57|1- ([H 

It is easy to see that entries of the Lax matrix (j47|) have the form (|55j) . So, we can use quadratic- linear algebra 
(|63[) in order to get bi-hamiltonian description of the restricted multiple three wave interaction system. The first part 
of the brackets between variables u,u*, bi,b* and Cj,c| may be directly restored from the brackets ([6"Tj) - ([6"2"|) . The 
remaining part has to be obtained from the compatibility conditions (|49|) . 

As a result the non-zero brackets with variables u and u* look like 



{u, = iu, {u*,bj}i = -%bjp 3 (\ + fi = ej), {u*,b*}i = -b*p 3 (X + p = ej) 



{u*,Cj}x = |Cjp 3 (A + /x = ej), {u*,c*}i = -- C jp 3 (\ + n = ej). 



The local brackets at j = 1, . . . , n are equal to 



{bj,Cj}i = 



2u'' 



{b . tC * h = l iik 



{bj,b*}i = 



i(bjC* — 2uej 



2u 



The non-local brackets at i ^ j read as 



ibiC *j ru l ib * h 3 r *i i6lC i r * mo ic i c J 



Other brackets are equal to zero, for instance 

{6 i ,6 J }i = {6 i ,c*} = {c|,c*} 1 = 0. 

Using r-matrix algebras (|48[) and it is easy to prove that integrals of motion Ij E { Jo, J%, ■ ■ ■ , J n , Ki, ■ ■ ■ , K n } 
(fT5|) are in the bi- involution with respect to the brackets {., .}o and {.,.}].: 

{I J ,I k } = {1^1^=0. 

The one of the main problems is that the main characteristic of the model, such as equation of motion, the Lax 
matrices and integrals of motion are invariant with respect to conjugation 

u <-> u* , bj <-> 6*, Cj <-> c*, (65) 

but the second brackets {., .}i do not invariant with respect to this transformation. 

The second problem is that the second brackets {., .}i are rational brackets and we have an obvious problem with 
quantization of these brackets. 

According to 0] there are different bi-hamiltonian structures for a given integrable system. So, we could try to find 
another bi-hamiltonian description of the restricted multiple three wave interaction system, which allows as to avoid 
these problems. 
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